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Survival Probability of the Neel State in Clean and Disordered Systems: 
an Overview 
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Abstract In this work we provide an overview of our re¬ 
cent results about the quench dynamics of one-dimensional 
many-body quantum systems described by spin-1/2 models. 
To illustrate those general results, here we employ a partic¬ 
ular and experimentally accessible initial state, namely the 
Neel state. Both cases are considered: clean chains without 
any disorder and disordered systems with static random on¬ 
site magnetic fields. The quantity used for the analysis is the 
probability for finding the initial state later in time, the so- 
called survival probability. At short times, the survival prob¬ 
ability may decay faster than exponentially, Gaussian behav¬ 
iors and even the limit established by the energy-time un¬ 
certainty relation are displayed. The dynamics at long times 
slows down significantly and shows a powerlaw behavior. 
For both scenarios, we provide analytical expressions that 
agree very well with our numerical results. 
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1 Introduction 

This overview describes our recent numerical and analytical 
results for the dynamics of quantum systems where exter¬ 
nal interactions with an environment are neglected mm 
I4ll5ll6ll7l . The focus is on the effects of the internal interac¬ 
tions and on the interplay between interaction and disorder. 
The system is initially in a non-stationary state very far from 
equilibrium. Our main goal is to understand what charac¬ 
terizes the dynamics of such isolated many-body quantum 
system. 

Nonequilibrium quantum physics is a subject that per¬ 
meates various fields of physics and chemistry, such as sta¬ 
tistical mechanics, condensed matter physics, molecular dy¬ 
namics, quantum information, and cosmology. It is also inti¬ 
mately associated with the development of methods to con¬ 
trol the dynamics of quantum systems, aiming at slowing it 
down or accelerating it. The subject is, however, much less 
understood than equilibrium quantum physics. 

We have been continually emphasizing in our previous 
works that the dynamics of a quantum system mi2i[^i5ii^ 
13 and also the new equilibrium reached by it ITTUSI depend 
not only on its initial state or on its Hamiltonian but on both. 
The effects of the two are intertwined. The system evolution 
depends on where the energy of the initial state falls with 
respect to the spectrum of the Hamiltonian and on how much 
spread out this state is in the energy eigenbasis. Both aspects 
depend on the details of the Hamiltonian, such as its density 
of states and the presence or absence of disorder. 

To connect our analysis with experimental studies, it is 
convenient to think about the dynamics in terms of “sud¬ 
den quenches”. The system is initially in an eigenstate of 
an initial Hamiltonian Hi. It is then taken far from equilib¬ 
rium by a sudden perturbation (quench) that changes Hi to 
a new final Hamiltonian Hp, initiating the evolution. There 
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is significant experimental interest in the quench dynamics 
of isolated many-body quantum systems. In nuclear mag¬ 
netic resonance (NMR), for example, the initial state can be 
prepared by applying a magnetic field to the sample, which 
once turned off starts the dynamics M- In this context, a 
system of particular interest for us is the crystal of fluorap- 
atite studied in solid state NMR MlOlll111 . The arrangement 
in this crystal is such that it can be treated for some time as 
a one-dimensional system of spins-1/2, which is similar to 
the systems considered in this work. Spin-1/2 models on a 
lattice can also be studied with trapped ions 012II13I and op¬ 
tical lattices 014lll5lfT^ . The latter offer several advantages, 
including high controllability, quasi-isolation, and flexibility 
in the preparation of the initial state. 

A state that has received much attention in experiments 
with optical lattices is the Neel state, partially because of 
its importance in studies about magnetism. The state is such 
that the polarization of the spin on each site alternates along 
a chosen direction. We use this state for our illustrations be¬ 
low. 

To quantify how fast the initial state l'f'(O)) = |ini) changes 
in time, we calculate the survival probability. 


F(f) = |(ini|e-'^^'|ini)|2 


a 


( 1 ) 


where C™ = (ini| i/a) is the overlap between the initial state 
and the eigenstates | i/a) of the final Hamiltonian and Ea are 
the corresponding eigenvalues. The survival probability has 
received several different names, such as fidelity, non-decay 
probability, and return probability. It gives the probability 
for finding the initial state at time t. Notice that if we know 
the envelope of the energy distribution of the initial state 
weighted by the components |C™p, that is. 


Pi„,(£)=LlC?5(£-£„), (2) 

a 

then by doing a Fourier transform we are able to obtain an 
analytical expression for F{t). This distribution is often re¬ 
ferred to as local density of states (LDOS) or strength func¬ 
tion. We use the first designation, but the reader should not 
confuse LDOS with the density of states. The latter is sim¬ 
ply the distribution of all eigenvalues of the Hamiltonian, 
without any reference to a specific initial state. 

When dealing with unstable systems, such as unstable 
nuclei, the decay is exponential, as observed experimentally. 
However, deviations exist. They are associated with the fol¬ 
lowing scenarios. 


1.1 Short Times 

By expanding Eq.([T]i, one sees that the initial decay has to 
be quadratic in time. This is the region of the quantum Zeno 
effect. 

F(f ^ 0) 1 - (3) 

where 

qni = ,/L|CrP(£a-£M)2 

y a 

= / L (4) 

y n/ini 

is the uncertainty in energy of the initial state and 

£,„i = (ini|i/f|ini)=£|c|,“|2£„ (5) 

a 

is the energy of the initial state with respect to the final 
Hamiltonian. Note that \n) denotes the basis vectors used 
to write the Hamiltonian Hf, the initial state |ini) being one 
of them. 


1.2 Intermediate Times: Strong Perturbation 

An exponential decay of the survival probability implies a 
Lorentzian LDOS. This is the regime of the Fermi golden 
rule. However, if the perturbation that takes Hi into Hf is 
sufficiently strong, being beyond perturbation theory, the de¬ 
cay of the survival probability can be faster than exponen¬ 
tial Il2imi4ll5ll^l71 . 

In realistic systems with two-body interactions, as the 
ones treated here, the fastest decay for a unimodal LDOS 
is Gaussian. In this case, the decay rate coincides with the 
width Oini of the LDOS: F{t) = exp(—Gaussian de¬ 
cays that continue beyond very short times were predicted 
for two-body random matrices (||23l and references therein). 
The novelty of our studies is the verification that this Gaus¬ 
sian behavior can indeed emerge in realistic systems and that 
it can persist until F(f) touches the saturation point. In addi¬ 
tion, we have shown that for some classes of initial states, it 
is straightforward to find Oini analytically i2im . 


1.3 Intermediate Times: Bimodal LDOS 

Another realistic case where the decay of the survival proba¬ 
bility can be faster than exponential occurs when the LDOS 
is bimodal. In this scenario, the decay at short times can 
reach the fastest velocity determined by the energy-time un¬ 
certainty relation, while the behavior at later times depends 
on the shape of each peak lH. 
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1.4 Long Times and Disorder 

Even if the system decays into the continuum, there is al¬ 
ways a lower bound, Ecut, in the spectrum. Taking this bound 
into account, the Fourier transform of the LDOS necessarily 
leads to a decay slower than exponential at long times El 
M- There has been analytical 019ll2OII2fl and experimen¬ 
tal II 22 I studies showing that the decay at long times should 
become powerlaw. The relation between the bounded en¬ 
ergy spectrum and a powerlaw fidelity decay at long times 
is shown below for two simple cases. 


7.4.7 Lorentzian LDOS 


Let us write the survival probability as F{t) = |A(f) f. In the 
absence of a lower bound, when Pini(7i) is a Lorentzian of 
width 77ni, the Fourier transform of the LDOS leads to the 
exponential decay of the survival probability, exp(—Tinif). 
We solve the integral, 

1 rini 


A{t) = 


f 


-iEt 


dE, 


Loo In (£i„i-£)2 + j-y4 
with residues. Since f > 0, we close the contour clockwise 
on the lower plane. 


^ r 

2n J-c 


-iEt 


r 

ini 


-E)+ iEiJl] [(£;„; - E) - 

^-iEt 

W^, 


dE 


-'Lm/2 


{E — Tiini) + i77ni/2] 

There is a pole at £ = Ei^i — i77ni/2. Taking into account 
the negative sign to 2ni, we obtain 


^(-27r0- 
2n^ ’ 




-iE,„,t , 


-/Tin' 


(6) 


On the other hand, when the lower bound does exist, 
A(f) can again be solved by replacing the integral with a 
contour integral in the complex plane E] , but the contour is 
now as follows. Choosing for convenience the energy bound 
to be Ecut = 0, the contour has the positive real energy axis, 
the arc of infinite radius running clockwise from the posi¬ 
tive real energy axis to the negative imaginary axis, and the 
negative imaginary axis from —i°° to the origin, that is. 


Pini {E ) £ 


-iEt 


dE = 

noo 

/ pini(£)e-‘®4£ 

Jo 

+ 

[ pM{E)e-‘^‘dE 


J arc 


/■O 

+ 

/ Pmi{E)e^ dE 

J —ioo 


(7) 

In general, the contribution from the integration along the 
circular arc vanishes. Thus, using E = —ie, we can write: 


Mt) = [ pii 
Jo 


i{E)e-‘^LE 


= j>Pmi{E)e 
= j) PmiiE)e 



‘^‘dE-if pini(-/e)e ‘^'de 

Jo 


= Ai(f) -|-A2(f). 


(8) 


The first term, Ai(f), depends on the poles in the fourth 
quadrant. As we saw in Eq. (|6l), it gives an exponential de¬ 
cay. It is the second term. 


A2(f) = -/ 


i-jL r 
2n Jo 


de ■ 


(£.n.-f/er+7I„V4 


that leads to the powerlaw behavior at long times, because as 
t becomes large, only small £ survives, so we can set £ = 0 
at the denominator. 


A2(/) —i 


r 

2n Jo 


de ■ 


1 

• cx: — 


+ Lni/^ f 


which implies that 

£(f ^oo)oc 1. 


(9) 


1.4.2 Gaussian LDOS 


Another simple example is the Gaussian distribution, in which 
case 

A{t) = ^ f (10) 

At very long times, the first exponential oscillates very 
fast, unless E is very small. Similar to what we did in the 
Lorentzian case, we can then set £ = 0 for the second expo¬ 
nential, 

A{t)'^ —oc(H) 

which again implies £(f —>■ oo) 

We note, however, that in our numerics for clean and 
disordered spin systems, the exponent of the powerlaw de¬ 
cay that we observe is < 1. Therefore, for our systems, the 
main origin of the algebraic decay of £(f) at long times is 
not the lower bound in the spectrum. Instead, this behavior 
is related to correlations between eigenstates. This is further 
discussed in Section 4, in the context of a disordered spin 
chain |[T]. 

This work is divided as follows. After describing the 
spin-1/2 models and the initial state considered in Section|2l 
we proceed to the two main parts of the paper. In Section[2 
we present situations where the decay of the survival prob¬ 
ability is faster than exponential. In Section 01 we show a 
powerlaw decay that emerges at long times in disordered 
spin chains. Concluding remarks are provided in Section|5] 
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2 Spin-1/2 Models and Initial States 

We consider a one-dimensional lattice of interacting spins- 
1/2 with an even number L of sites and two-body interac¬ 
tions. The Hamiltonian is given by 

H='^hkSl + 
k=l 

k ^ ' 

+A (SISU^ + ^k^k+2 + ^SIS1^2) ■ (12) 

In the equation above, we consider the Planck constant h = 
1, are spin operators, /zj. are the Zeeman splittings of 
each site k,J=l sets the energy scale, and A is the anisotropy 
parameter. The sums in the second and third lines run from 
k = 1 to L for periodic boundary conditions and from k = 1 
to L — 1 for open boundary conditions. The total spin in the 
z-direction, = Y^k^l, is conserved. We work with the 
largest subspace, = 0, of dimension ^ = L!/(L/2) !^. 

For the clean system investigated here, the magnetic field 
along the z axis is constant, /z*. = /z. The disordered system 
is characterized by random static magnetic fields, where the 
amplitudes /z*. are random numbers from a uniform distribu¬ 
tion [—/z,/z]. We also look at the case where a single defect 
exists, that is, only one site has /z*. different from the others. 

In the absence of disorder (/z^. = /z) and of next-nearest- 
neighbor (NNN) couplings (A = 0), the Hamiltonian is solv¬ 
able with the Bethe ansatz ll24ll . This clean Hamiltonian with 
only nearest-neighbor (NN) couplings is referred to as XXZ 
model. Disorder 025II26I|Z71 . even a single defect 02811291 . 
or NNN couplings 029113011311 can take the system into the 
chaotic domain. 

Initial State. We take as initial state one of the so-called site- 
basis vectors (also known as computational basis vectors). 
They correspond to states that are on-site localized, so the 
spin on each site either points up in the z-direction or down. 
The initial state selected is the Neel state, which is given by 

iNs) = intinnn...). (b) 

The scenario where a site-basis vector evolves accord¬ 
ing to the Hamiltonian given by Eq. (fT2t is equivalent to a 
quench where the initial Hamiltonian is the Ising part of the 
total Hamiltonian, Hj = Y.,k^k^k+v Hamilto¬ 

nian is Hf ~ H from Eq. (fT2l . The perturbation to change 
the Ising Hj to Hf, where couplings in the .ry-plane also ex¬ 
ist, is very strong, being beyond perturbation theory. As a 
consequence, the initial decay of the survival probability is 
expected to be faster than exponential. 


When the initial state is a site-basis vector, it is straight¬ 
forward to obtain its energy [Eq. (|5]l] and its energy un¬ 
certainty Oijji [Eq. dUi]. As indicated by those equations, the 
energy is simply the diagonal element of the Hamiltonian 
matrix written in the site-basis and is the sum of the 
squares of the off-diagonal elements in the row of that cho¬ 
sen initial state. To obtain ■, we just need to count how 
many site-basis vectors are directly coupled to the chosen 
initial state. 

3 Clean Systems: Faster than Exponential Decays 

We analyze two examples where the decay of the survival 
probability is faster than exponential: when the LDOS is a 
single Gaussian and when it is composed of two well sepa¬ 
rated Gaussians. 

3.1 Gaussian Decay 

Let us analyze first the case where the final Hamiltonian de¬ 
scribes an anisotropic clean system with open boundary con¬ 
ditions, NN and NNN couplings: = /z, A = 1/2, A = 1. 
The density of states, as shown in Eig.[T](a), has a Gaussian 
shape. This is typical of systems with two-body interactions 
II32II33I . although the distribution is not necessarily symmet¬ 
ric 03411351 . The Gaussian form implies that most states con¬ 
centrate in the middle of the spectrum. This is the region 
where strong mixing can occur and where the eigenstates 
are therefore more delocalized. As a consequence, an initial 
state with energy Q close to the center of the spectrum 
decays faster than a state with closer to the edges. 

The initial state considered is the Neel state [Eq. ([T3]) 1. 
Its LDOS is also Gaussian as seen in Eig.[T](b). This is an in¬ 
dication that this initial state has reached its maximum pos¬ 
sible spreading, which is limited by the density of states. 
Such maximum spreading is due to the strong perturbation 
that takes Hi to Hf. The LDOS is centered at the energy 
obtained from Eq. which for the Neel state is 

^open,clean ^ 1) + (L-2)A]. (14) 

The width of the LDOS corresponds to the uncertainty in the 
energy of the initial state [Eq. (|4]i], which for our initial state 
is 

^open,clean^A^/^ (15) 

Eor the chosen parameters of the final Hamiltonian, 
is close to the middle of the spectrum, which explains why 
the LDOS is well filled. Eor states with energy closer to 
the edges of the spectrum, the LDOS becomes more sparse. 
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Fig. 1 Density of states (a) and (d); LDOS (b) and (e); and survival probability decay (c) and (f). Open boundary conditions; L = 16. Clean 
Hamiltonian with NN and NNN couplings, 4 = 1/2, A = 1 (a,b,c) and the same Hamiltonian with an additional defect e = 15 on site L/2 (d,e,f). 
Initial state is the Neel state on (b) and (c) and an equally weighted superposition of two Neel states on (e) and (f). The solid lines in (b) and (e) 
are the Gaussian distributions with width (S). In (c) and (f), the solid lines are the numerical results, while the circles coiTespond to the analytical 
expressions, Eq. GU and Eq. Go), respectively. Dashed lines give the saturation points [Eq. E)]- 


Some examples of this latter case may be found in ||2ll3ll4l 

Ea. 

The decay of the survival probability is shown in Fig. [T] 
(c). The behavior is Gaussian (solid line) all the way to the 
saturation point (dashed line) and it agrees very well with 
the analytical expression (circles), 

^’(0=exp(-(j2jif^), (16) 

with CTini given by Eq. (O. 

The decay of the survival probability eventually satu¬ 
rates because we deal with finite systems. In a system with¬ 
out too many degeneracies, the off-diagonal terms in the sec¬ 
ond term on the left-hand side of the equation 

F{t) = Y.\cT\'+ L 

a 

averages out, leading to the infinite time average, 

^ = LlC‘l" = IPRmi- (17) 

a 

One sees that the saturation point depends only on how 
much spread out the initial state is in the energy eigenbasis. 
If many eigenstates | xj/a) contribute to the evolution of |ini), 
then there are many very small components |C™|^ and the 
saturation point is low. 

In Eq. (fTTl i. IPR stands for inverse participation ratio. 
This is one of the most common quantities used to mea¬ 
sure the level of delocalization of an arbitrary state |^) = 
Y.jCj\(l)j) in a certain basis |0), 

IPR|0=E|C/l"- (18) 


When the state coincides with one of the basis vectors, then 
IPR= 1. The more delocalized the state is in the chosen basis, 
the smaller the value of IPR is. The minimum value for sys¬ 
tems with time reversal invariance is reached by the eigen¬ 
states of real and symmetric full random matrices. In this 
case, the states are random vectors, leading to IPR^ 3 

m- 


3.2 Cosine Square Decay: Lower Bound from the 
Energy-Time Uncertainty Relation 

We consider again a final Hamiltonian with NN and NNN 
couplings, open boundary conditions, 4 = 1/2, A = l, but 
now not all sites have /i<. = h. Site Ljl is the defect site with 
hi^l 2 = h-\-e, where e is an excess on-site energy. Eor e ^ 1, 
the defect effectively breaks the chain in two. The spectrum 
ends up having two sets of eigenvalues. One set corresponds 
to the states that do not have an excitation on the defect site, 
so they have low energies, and the other consists of the states 
with an excitation on the defect site, so they have large ener¬ 
gies. The resulting density of states is a bimodal distribution 
with two well separated Gaussians, as shown in Eig.[T](d). 

To study the effects of these two peaks on the dynamics, 
we choose as initial state a superposition of two Neel states, 
so that one has an excitation on the defect and the other does 
not, 

1,,,^, I • • •) +1 AUUUtit ■■■) 

P^^/sup — /- 


(19) 
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The LDOS for this superposition is shown in Fig.[T](e). 

As expected, it is also bimodal with one peak centered at 
^open,clean _ ^open,clean ^ 

peaks have the same widths given by KEa.fTsTlI. 

The decay of the survival probability, shown in Fig. [T] 
(f), agrees very well with the Fourier transform of the two 
Gaussian peaks. 


F{t) = cos^ 



/ open,clean\2^2 
«s n 


( 20 ) 


The decay at short time is dominated by the cosine square 
part of Fq. (l20l i. which later leads to revivals. The envelope 
of these subsequent damped oscillations is determined by 
the Gaussian part of Fq. (l20l i. 


Lower Bound Established by the Energy-Time Uncertainty 
Relation. The cosine square decay during t < n/e is the 
fastest decay of the survival probability allowed by the energy¬ 
time uncertainty relation. 

The lower bound for the 
I44ll45]l46ll . 


F(f) >cos2(ainif), 


( 21 ) 


can be derived from the Mandelstam-Tamm uncertainty re¬ 
lation. 


1 

> - 


d{A) 

dt 


The uncertainty in energy, Gh, for a non-stationary state co¬ 
incides with (Tinp If A is the projection operator on the initial 
state, A = |ini)(ini|, then (A) =F{t) and =F(t) — F(t)^. 
Thus 




dF 

dt 


To solve this equation we can follow Ref. Il43l and write 
F{t) = cos^ 0, so that 


Ojni COS 0 sin 0 > - 


dO 

—2 cos 0 sin 0 — 
dt 


dO 


dt 
arccos 


< (Jini =» |0(f)| < (Tinif 




■\Vm) 

Since the arccos is a strictly decreasing function. 


arccos (^s/F{t)'j > a-^nit ^ \jF{t) > cos((Jinif), 

from where Fq. dZB follows. Notice that this expression is 
valid only while the cosine is positive. In our case, where the 
FDOS is composed of two well separated peaks, ^ e/2. 


4 Disordered Systems: Powerlaw Decays 

For the relatively small clean systems studied above, it is 
not straightforward to identify the powerlaw decay at long 
times. In the case of Fig. [T](c), for instance, we need to do 
a time average of F{t) to notice that there is indeed an alge¬ 
braic decay for 2 <t < 100, before the the survival probably 
starts to simply fluctuate around its infinite time average. 
The powerlaw exponent of this decay is very close to 1, that 
is F{t) 

To explore the onset of smaller powerlaw exponents, we 
now add disorder to the system, which slows down its dy¬ 
namics. We consider the disordered XXZ model described 
by an anisotropic Hamiltonian (A = 0.48) with periodic bound¬ 
ary conditions, NN couplings only (X = 0), and random 
static magnetic fields, where /ij. are random numbers uni¬ 
formly distributed in [—h,h]. 

In a non-interacting system (A = 0), the presence of dis¬ 
order localizes the excitations. This is the scenario of the 
Anderson localization El , which has been extensively stud¬ 
ied O48ll49ll50l and also tested experimentally, more recently 
in two-dimensional ultracold gases with speckle disorder EH. 

It had been conjectured that localization should persist 
also when interactions would be taken into account II47II52II . 
This was confirmed with perturbative arguments in 053II54I 
and rigorously in ESI. So far, this so-called many-body lo¬ 
calization (MBF) has been tested in one experiment with 
cold atoms in optical lattices Il56l . Here, we study the evo¬ 
lution of the Neel state as the system approaches the MBF 
phase. 

When /r ~ 1, the disordered XXZ model is chaotic ll25l 
12^2711281 . so the eigenstates close to the middle of the spec¬ 
trum are very delocalized in the site-basis, being similar 
to random vectors. Their inverse participation ratio is in¬ 
versely proportional to the dimension of the symmetry sec¬ 
tor to which they belong. Analogously, initial states corre¬ 
sponding to site-basis vectors with energy close to the center 
of the spectrum of a chaotic Hp have IPRini ^ * ■ 

As h increases, the eigenstates jy/a) become more lo¬ 
calized, sampling only a portion of the Hilbert space, so 
IPRa where £>2 < 1. Equivalently, the initial site- 

basis states get less spread out and IPRini with 

£>2 < i- The exponents D 2 and D 2 are known as general¬ 
ized dimensions. For one-dimensional systems without in¬ 
teraction, they coincide, D 2 = £> 2 ^ as shown in Ref. E3. We 
verified that the same holds for our interacting system IH . 


Anderson Localization and Powerlaw Decay. In studies of 
the Anderson localization, it has been shown that the decay 
of the survival probability averaged over random realiza¬ 
tions, {F[t)), becomes powerlaw at the critical point, with 
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an exponent that coincides with D 2 058II59IIM1 . This gener¬ 
alized dimension is obtained from the scaling analysis of the 
inverse participation ratio IPRini- 

The subscript “2”, appearing in D 2 and £> 2 , is used to dis¬ 
tinguish from generalized dimensions obtained from scaling 
analysis of Y,a with q£^2, which are not considered 

in this work. They are important when investigating multi¬ 
fractal features of the states. 

The coincidence between the powerlaw exponent and D 2 
may be understood from the expression of the survival prob¬ 
ability as follows. Introducing the identity 

/ oo 

do}e'^'8{Ea-Efi-o}) (22) 

into {F{t)) fromEq. ([T]i, we obtain: 

(^( 0 ) = 

= J dcoe‘“'C{co), (23) 

where C{co) is the correlation function given by ||62l, 

C(®) = (E„,/3 l4lV“niP5(£„-£/i -0))^ . (24) 

The correlation function C(a)) quantifies the overlap between 
two eigenstates with energy difference [isa — E/i | ■ In the crit¬ 
ical regime, C{(o) decays as a powerlaw 15911631 : for small 
CO and large it scales as I59II60I . 

C(® ^0)oc ^0)^2-'. (25) 

At large f, due to the rapidly oscillating exponential term in 
Eq. (l23l l. the integral is dominated by small co. If we then 
substitute Eq. (|25]) into Eq. (|2^ . we obtain 

(26) 

We see that D 2 measures how much correlated the compo¬ 
nents of the initial state, are and equivalently how 

much correlated the eigenstates are. In a chaotic system, 
where the states are uncorrelated random vectors, D 2 ^ \ 
and C{co ^ 0) As the disorder increases, the states 

shrink, correlations build up, and £>2 becomes smaller than 
1 . 

We note, however, that the comparisons between the gen¬ 
eralized dimension and the powerlaw exponent have in fact 
been carried out for the time-averaged survival probability, 
defined as 

C(f) = I f{F{T))dT. (27) 

t Jo 

This is a way to smooth the curve of the survival probabil¬ 
ity. In order to reduce also the fluctuations in the values of 
IPRjni, the scaling analysis is often done with the so-called 
typical inverse participation ratio, IPR'^p = exp((lnIPRi„i)). 
The scaling analysis of IPR'^' gives Z)™- 


Many-body localization and powerlaw decay. Recently, we 
verified that the correspondence between the powerlaw ex¬ 
ponent and the generalized dimension holds also in the pres¬ 
ence of interactions [T] . The agreement is excellent for both 
C(f) and {F{t)) when the disorder is small and the system 
sizes are large. As the disorder increases, the largest sys¬ 
tem sizes available for exact diagonalization, L = 14,16, are 
still too small. In this case, oscillations are observed before 
the powerlaw behavior, with D 2 now capturing the decay of 
these oscillations. There is agreement between D '2 and the 
powerlaw exponent of C(f), but for a short time interval. 

In Ref. Q , we performed averages over realizations and 
also over initial states. The latter was a set of site-basis vec¬ 
tors with the closest energies to the middle of the spectrum. 
In the present work, the initial state is always the Neel state, 
so the averages are performed only over realizations. The 
averaged energy of the Neel state is 


This value is still far from the edge of the spectrum, but 
it is not at the center of the spectrum either, so it should be 



Fig. 2 LDOS for one disorder realization (a) and (c) and survival 
probability averaged over 10^ disorder realizations (b) and (d) for 
h= 1.5 (a,b) and h = 3 (c,d). Periodic boundary conditions; A = 0.48, 
A = 0, L = 16. Initial state is the Neel state. Solid lines in (a) and (c) 
are Gaussians centered at (Fjj'""' *”) [Eq. (2^ 1 with width given by 
Eg. 1291 . In (b) and (d), the bottom solid curve gives {F{t)) and the 
top solid curve coiresponds to C(f). The bottom dashed line gives 

_f^yp 

and top dashed line gives . Circles indicate the Gaussian decay: 
exp [—■ Thick horizontal lines represents the saturation 
points. 
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easier to localize this state than the states that we considered 
in m . Notice that it is the absence of NNN couplings that 
pushes the energy of the Neel state away from the middle of 
the spectrum [compare the equation above with Eq. ([I4li 1. 

InFigs.|2](a), (c), we show the LDOS for one disorder re¬ 
alization for two values of the disorder strength, !i = 1.5 and 
h = 3, respectively. The distribution becomes visibly more 
sparse as h increases, although the width does not change, 
since it is independent of h, 

/_closed, diso\ fT 

As a consequence, the short-time dynamics of the survival 
probability, which depends only on the width of the LDOS, 
is identical for both values of h, as shown in Figs.|2](b) and 
(d). 

The bottom solid curve in Figs. |2](b) and (d) indicates 
{F{t)) and the top solid one represents the average C(f). 
Both {F{t)), and C{t) reach the same saturation point and 
this values naturally increases as h increases. 

The initial Gaussian decay of {F{t)) is followed by damped 
oscillations and then finally saturation. It is interesting that 
at short times, around f ~ 2, {F{t)) reaches values signifi¬ 
cantly lower than the saturation point. This is a result that 
deserves further investigation. Site-basis states that, like the 
Neel state, overshoot their decay seem to be behind the emer¬ 
gence of time intervals where (T(f)) remain below the satu¬ 
ration point, as it observed in Fig. 1 of Ref. IT]. 

In Figs.|2](b) and (d), the algebraic decay (the bot¬ 
tom dashed line that is closest to {F{t))) matches the enve¬ 
lope of the decay of the oscillations of {F{t)). This occurs 
even for !i= 1 .5, which is in contrast with our results in m, 
where the oscillations were minor for small disorder and D 2 
agreed very well with the powerlaw exponent. The reason 
for this discrepancy is that even close to the chaotic regime 
of the disordered XXZ Hamiltonian, the Neel state is not as 
delocalized as the states we studied in m , because of its low 
energy, so oscillations appear already at small h. In fact, for 
this disordered system, the Neel state does not seem to be 
able to reach a diffusive behavior, where {F{t)) f^*, hav¬ 

ing instead D 2 ,D '2 < 1 for all values of li> 1 [see also Fig. [3] 

(d)]. 

The oscillations of {F{t)) are substituted by a powerlaw 
behavior when the average C(f) is considered [Figs. |2](b) 
and (d)]. Despite holding for a relatively short time interval, 
there is reasonable agreement between the numerical curve 
and t 2 (the latter corresponds to the top dashed line that 
is closest to C(f)). 

We expect the agreement between the exponent of the 
powerlaw decay and the generalized dimensions to improve 
and the powerlaw behavior to persist for longer times for 
system sizes larger than L = 16. Good indication for this 



Fig. 3 {Fit)) (a) and C{t) (b) for h— 1.5; L = 12 (top solid curve) and 
L = 16 (bottom solid curve); initial state is the Neel state. Dashed lines 
give (a) and (b). For (c): Scaling analysis of InIPR* vs 

ln,/T for h= 1.5 (circles) and h = 3 (squares). Error bars are standard 
deviations over 10^ values of lnIPRi„i. For (d): Dj’’ vs disorder strength. 
Solid line is a fitting curve. All panels: periodic boundary conditions; 
A = 0.48, A = 0. 

claim is found in Fig. [2(a) and (b), where we compare the 
survival probability for L = 12 and L = 16 for h= 1.5. Fig- 
ure|3](a) shows {F{t)) and Fig.[3](b) gives C(f). It is evident 
that as the system size increases, the duration of the oscil¬ 
lations of {F{t)) and of the algebraic decay of C(f) stretch 
out. One can see that for L = 16, the damping of the oscil¬ 
lations coincides better with than for L = 12. Similarly 
matches C{t) when L = 16, but this is not the case for 
L= 12. 

In Fig. 12(c), we show two examples of the scaling anal¬ 
ysis that we perform for IPR-^'j to obtain The slope nat¬ 
urally decreases as the disorder strength increases, while the 
standard deviations increase with it. The value of the gen¬ 
eralized dimension becomes less precise as the system ap¬ 
proaches the MBL phase. 

The dependence of D* on the disorder strength is de¬ 
picted in Fig. 12(d). The decay clearly slows down for h > 
5, but since the system sizes considered are small and fi¬ 
nite size effects become more relevant for larger disorder 
strengths, we avoid speculating on this behavior. 


5 Concluding Remarks 

We studied the survival probability of one-dimensional sys¬ 
tems of interacting spins-1/2 in the absence and presence of 
disorder. The initial state considered was the Neel state. We 
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showed that its short-time dynamics is characterized by the 
shape and width of the LDOS, and it does not depend on the 
strength of the disorder. The long-time dynamics depends on 
how well filled the envelope of the LDOS is. As the disor¬ 
der strength increases, the LDOS becomes more sparse and 
the decay at long times slows down significantly. The decay 
becomes powerlaw with an exponent that coincides with the 
generalized dimension. 

We have numerical results for various observables of ex¬ 
perimental interest, such as magnetization and spin-spin cor¬ 
relations. In the presence of disorder, they also show alge¬ 
braic decays. A natural extension to the present work, which 
we are studying, is the search for analytical expressions also 
for these observables. One expects the expressions obtained 
for the survival probability to help in this direction. 
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